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S ' Abstract 

Every polyhedral cone can be described either by its facets or by its extreme rays. 
Computation of one description from the other is a problem that can be very complex, 
i.e. one encounter the combinatorial explosion. We present here several methods which 
allow us to use symmetries for the computation of extreme rays and facets of cones. I 
give an algorithm for computing adjacencies using a linear programming approach. I give 
a complete description of the method of incidence, adjacency decomposition method and 
• of the subcone method used for computing with symmetric cones. 

d ■ 

1 Introduction and motivation 

>- [ A cone in is a set stable by addition and multiciplication by A G M+. A polyhedral cone is a 
2 ■ cone defined by a finite number of linear inequalities. It is a theorem of Minkovski that those 
polyhedral cones can also be described as positive span of a finite number of rays. A ray is 
^ . called an extreme ray if it is not the sum of two different rays. An inequality is called a facet 

I if it cannot be expressed as a sum of two different inequalities. 
O I So, there is essentially two description for a cone, either by extreme rays or by facets. The 

^ • main problem of this paper is given a description by facets to obtain a description by extreme 
d • rays. The dual problem of computing a description by facets from a description by extreme 
S . rays is equivalent to the previous one. 

j> ' In many mathematical, physical, and practical problems we have a description by facets (or 

' extreme rays), the interesting objects being the extreme rays (or facets). So, the problem of 
translating the descriptions is a very interesting one and many programs exist that do it (cdd 
^ ([Fu95]), porta \ Irs ([Av]), etc). 

But, in many cases one encountered the so called, combinatorial explosion, i.e. while the 
description by facets (or extreme rays) remain simple, the dual description has a huge number 
of extreme rays (or facets). 

These cones have in many case some symmetry groups. Our goal in this paper is to present 
several methods for computing extreme rays and facets of cones using their symmetry group. 
Two classical algorithms exist for computing extreme rays: 

1. The Double description method, also called Fourier Motzkin method, consists of taking 
an initial list of n facets of the cone and computing extreme rays of it (it is a simplicial 
cone so this is immediate), then to add the remaining facets one by one in order to have 
the list of extreme rays. The method is called Double Description since at any time of the 
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process we know the description by facets and by extreme rays. This method is the one 
of cdd by Komei Fukuda and Avis (see [FM] and [FP]). This method has some variants 
(see for example [VB]) 

2. The Pivoting method consists of taking an extreme ray of the cone and finding extreme 
rays adjacent to it. This method is generally inefficient for extreme rays with high de- 
generacy that occur in combinatorial optimization. 

The cones and polytopes appear in many context like Metric Geometry (metric cone, quasi- 
metric cone, super-metric cone (see [DeLa97], [DDOl], [La], [Gr92]), Max-cut problem (cut 
cone and cut polytope, see [DeLa97] and [D-L]), Geometry of numbers (Voronoi polytope and 
Voronoi Domain, see [Jaq]), Combinatorial optimization (Traveling Salesman Polytope and 
Linear Ordering Polytope, see [ChRe96]), Game theory (solitaire cones, see [AD]), Physics (see 
[BCLM] and [Er]). 

There are two main kinds of symmetry groups: Permutation groups like Sym{n) arising in 
combinatorial optimization and matrix groups arising in geometry of numbers (see [Jaq]). 

Note also that the problems of finding the size of an orbit and of finding canonical repre- 
sentatives may be very difficult (see for example [LN], [Jaq]). Althought we use the simpler 
language of cones here, the methods apply also to polytopes. 

2 General infos on cones and polynomially solvable ques- 
tions 

A cone in M" can be defined in several ways, first as an intersection of a finite number of half 
spaces: 

n 

X eW such that ^ ai^kXi > 0, VI < < m 
1=1 

or as a set spanned by vectors: 

p 

X eW such that 3Xi > with a; = ^ XiVi . 

1=1 

We assume that the cones have full dimension n. Given v e K", the inequality J2i=i ''^i^i ^ 
is said to be valid for C, if it holds for all x e C. Then the set {x e C\ Yl^=i "^i^i — 0} is called 
the face of C , induced by the valid inequality Y17=i "^i-^i > 0. A face of dimension n — 1 is called 
a facet of C; a face of dimension 1 is called an extreme ray of C. 

Two extreme rays of C are said to be adjacent, if they span a two-dimensional face of C. 
Two facets of C are said to be adjacent, if their intersection has dimension n — 2. The 1-skeleton 
graph of C is the graph whose nodes are the extreme rays of C and whose edges are the pairs 
of adjacent extreme rays. The n — 1-skeleton graph of C is the graph with node set being the 
set of facets of C with an edge between two facets if they are adjacent in C. 

Some face conditions may be useless to the cone description by face, they are called redun- 
dant. Some vector may also be redundant. In fact, we can find a finite subset of non- redundant 
faces of the cone C called facets. The non-redundant rays are called eodreme rays. 

Given a cone C we can define C* as the cone whose extreme rays are facets vector of C, the 
facets of C* will be the extreme rays of C. So, every method or theorem given here as a dual 
version with extreme rays and facets exchanged. 



The set of faces of a cone C is a lattice denoted by L{C) under the inclusion. Given a set JF 
of faces a face named sup JF is defined as the smallest face containing the elements of JF, while 
inf is defined as the biggest face contained in every element of J^. 

The k skeleton of C is the set of k dimensional faces of C. Two /c-faccs Fi, F2 are adjacent if 
sup {Fi, F2} has dimension k + 1 and inf {Fi, F2} has dimension k — 1. The k-skeleton graph is 
the graph with node set being the set of A;-faces of C with an edge between two /c-faces if they 
are adjacent in C. 

The linear programming way of testing adjacencies 

1. First take a set ^ = {fi}iei of inequalities fi{x) > and another inequality /, the 
question is: Is that / redundant? reply is to consider the linear programming problem: 

minimize f{x) in the cone defined by J- 

if the minimum is then / is redundant (i.e. it can be expressed as a sum with non- 
negative coefficients of elements of J-). If the minimum is —00 then / is is non- redundant. 

2. Claim Given a non-redundant set JF of facets of a cone C and (/, g) a pair of facets then 
{f,g) are non-adjacent if and only if there exists Aj > such that 

Kfi = A// + Xgg 

with Xf > 0, Ag > 0. 

3. Using the above claim we have {f,g) non-adjacent if and only if 

infxGC f{.x) = 
with C = {x such that u{x) > 0, u e (JF - {f,g}) U {-g)} 

This method allow us to compute adjacencies of facets of a cone defined by facets without 
knowing its extreme rays. It is so very efficient. 

The linear programming method for computing k-faces 

1. We assume here that a cone is defined by a set STZ of extreme rays. 

2. The linear programming method allow us to compute edges from the set of extreme rays, 
i.e. we find from the set of extreme rays the set of edges. Prom this set of edge one would 
hke to compute the set of adjacent edges, i.e. we want to compute the set of edges that 
span a 3 face of the cone C. 

More generally generally the computation of adjacencies of fc-faces amount to the com- 
putation of elements of the k + 1-skeleton. 

3. Note that the edges are 1-simplex but the clement of the elements of the 3 skeleton can 
be nonsimplicial. This fact can be seen in several ways: 

(a) The 3 skeleton of icosahedron is made of triangles while for dodecaedron it is made 
of pentagon. 

(b) n = 3 is the first dimension where cones can have an arbitrary high number of 

extreme rays. 

(c) In a list of 3 rays that are linearly dependent, at least one of them is redundant, 
while there exist list of 4 rays that are linearly dependent but all irredundant. 



4. Linear Programming way of testing the adjacency: Let Fi and F2 be two A;-faces, we 
assume that inf{Fi, F2} has dimension k—1 (i.e. it is one of the /c — 1 faces found at the 
last step). The process is as follow: 

(a) We determine the set £ of extreme rays e such that e G Vect (FiU, F2) 

(b) sup{Fi U F2} has dimension superior to A; + 1 if and only if there exist Ag' > and 
Ae > such that 

K'e' = Aec 

e'esn-s ee£ 
this correspond to the linear programming 

infa;£c es.x = 

with C — {x such that e.x > Ve e £TZ — £ and e.x < Ve e £ — {es}} 

3 The method of Symmetry 

Consider a cone C with a linear symmetry group G acting on it. If F is a facet then the 
symmetry group of the facet is 

Symp = {g E G such that g{F) = F} 

is a subgroup of G satisfying to ISympl x \0f\ = 1^*1 with Op the orbit of F. Furthermore if F 
and F' belong to the same orbit then Syrup and Syrup' are conjuguate subgroups of G. The 
conjugacy classes of Sym{n) are 



n 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Conjugacy classes 


3 


5 


7 


11 


15 


22 


30 


42 


56 



If if is a subgroup of G then we define Ch = {d & C,\/h G H, h{d) = d}. The cone Ch has a 
dimension smaller than n. If C is defined by extreme rays then the extreme rays of Ch can be 
found easily. 

Theorem 1 If F is a facet of C having H as a symmetry group then F is a facet of Ch- 

We setUp = {dEC s.t. F{d) — 0}. If F has H symmetry, then Up has also H as a, symmetry 
group. 

Now define for any ray of C 

SU{d) = h{d) 

heH 

Denote SU{X) = {SU(d) s.t. d e X} for any X C OMCUTa. Then SU{Up) is a set 
of symmetric elements, which are incident to F. Moreover, since Ch = SU{C), we have 
SU{Up) C C. Wc denote n the dimension of C and m the dimension of Ch- By hypothesis F is 
a facet, so Up has dimension n — 1. The mapping d — SU{d) decrease dimension by at most 
n — m; so, we get that SU{Up) has dimension m — 1, i.e. F is a facet of Ch- 



Nevertheless the facets of the cone Ch are not necessarily facets of C, i.e. after computing 
facets of Ch, we need to check if they are indeed facets of C. 

Using this method we can found the facets of a given cone having a given symmetry. Since 
the number of conjugacy classes of a group is not large we can compute orbits having a given 
non-trivial symmetry. 

The cone CUT-j has the property that all its facets have a non-trivial symmetry group but 
this is a particular case and in general we can't expect to find all facets in this way. 



4 The method of incidence 



Let C be a cone of dimension n with a symmetry group G and a set !F of facets partitionned 
into orbits {Oi, . . . , Op} with representatives {Fi, . . . , F2} 

If e is an extreme ray then it is incident to at least one facet of C. By the group action, 
we can assume it is incident to a representative of orbit O^. This incidence yield a linear 
equality and this equality allow us to find a cone of dimension n — 1 of which e is an extreme 
ray. This idea can be extended to faces of dimension n — k with k small. Each such faces is a 
cone of dimension n — k. 

This idea was used by Grishukhin in [Gr92] on METj and by Christof and Reinelt in 
[ChRe96] and by [DMP] for the metric cone METg. The method is the following: We find a 
set of orbits of n — A; faces of C using the linear programming method of the preceding section. 
Then for each of these n — k-iaces we find their extreme rays using classical program like cdd 
([Fu95]). 

One problem of the method is that a given extreme ray belong to many n — /c-faces so the 
orbit of extreme rays found will belong to different orbit of n — /c-faces. 

5 Adjacency decomposition method 

Take a cone C of dimension n defined by a set £ of extreme rays of C. Now consider the set 
Sp of extreme rays incident to a given facet F. The facet F can be considered as a cone Cp of 
dimension n — 1 generated by the set Sp of extreme rays. 
Furthermore the mapping 

(L(C),C) ^ {L{Cp),(l) 
G ^ GnF 

is a lattice morphism. In particular the facets of the cone Cp correspond to the ridges included 
in F. It is computationally easy to map a facet of Cp to a ridge R of C, to find the facet G of 
C such that R = G H F and so, to build facets of C. 

So, we can from a facet of C find the adjacent facets to F, the complexity of the computation 
being the one of finding facets of the cone Cp defined by its extreme rays, this complexity is 
related to the incidence number of F. If it is equal to n — 1 then the facet is simplicial and this 
is trivial, if the facet has incidence larger than n — 1 then this becomes more difficult and we 
encounter again the combinatorial explosion while on a lesser form. 

The algorithm is so: 

1. Take a list of orbits of facets of the cone Cp, find their incidence numbers and orbit size. 

2. Take a facet from the list which is non-treated and find adjacent facets to it as indicated 
above 

3. Extract from this list of facets the orbits and update the hst of orbits. 

4. Go to step 2 if there remain some non-treated facets. 

This method was applied by Deza, Fukuda, Pasechnik, and Sato in [DFPS] for the metric 
polytope METs and METg, by Ghristof and Reineh, in [ChRe96] for CUTs, CUTg, the hnear 
ordering polytope P|q and the Traveling Salesman Problem P^sp and P^sp, by Deza, Dutour 
and Pantaleeva in [DDP] for OMCUT5, OMCUTq and QMETq, by Deza and Dutour in [DDOl] 
for supermetric cones SMET^'^ . 

The result of this method fall in three cases; 



1. The method find all orbits of facets, this is the case of OMCUT^, P^sp ^ot example 

2. The method find a set of orbits of facets which is conjectured to be complete, this is the 
case of HCUT^, HCUT^, CUTs, MET^, SMET^''^ 

3. The method give only an estimate of the number of orbits, this is the case of METg, 
CUTg, P^sp: OMCUTe, SCUTq'^. 

So, this method allow us to draw a frontier of the combinatorial explosion. In the second 
case, we find facets adjacent to orbits except for a few ones that have a large number of incident 
extreme rays, example 

1. For CC/Ts, the triangle facet, the pentagonal facet and the seven gonal facets have a too 
large incidence number (see [ChRe96]) 

2. For HCUTy and HCUTq, the non-negativity facets and the simplex facets. 

3. For METs the cuts (see [DFPS]). 

4. For the Travehng Salesman Polytope, the Linear Ordering Polytope see [ChRe96] 

The complexity of these facets is usually extremely high, but one may want to prove that the 
found description are complete. 

The mathematical object that appear is the quotient graph of the skeleton graph with 
elements, the orbits of facets. Two orbits Oj and Oj being adjacent if there is a facet 
Fj e Oi and a facet Fj e Oj such that Fj and Fj are adjacent facet in the cone C. 

The mathematical conjecture is that the quotient graph with those orbits removed is still 
connected. As far as we know no result of this kind as ever been proved. The only result 
available is the Balinski theorem which asserts that the skeleton graph of a rf-polytope is d 
connected, i.e. the removing of d elements leave it connected. This result could not be applied 
in our cases since the size of the problematic orbits is larger than the dimension of the underlying 
space. 

Another approach is possible: These orbits of facets have a very simple expression and 
a large symmetry group (for example the symmetry group of the triangle facet of CUTn is 
Z2 X Sym{n — 3)) so, the preceding approach can be applied recursively to the cone Cp- This 
approach succeeds with HCUT^, we intend to apply it to the cut cone CUTg. 

This approach is the one of Jaquet in [Jaq]. [NEED TO CHECK THIS, THE JAQUET 
METHOD SEEMS QUITE INVOLVED AND SMART, HE REWRITES EVERYTHING] 

In order to run the adjacency decomposition method, one needs to have an initial set of 
facets, among the possible method are: 

1. Use of linear programming to generate vertices 

2. Special kind of facets known by the theory of particular cones (zero-extension for metric 
cones for examples) 

3. Partial cdd output obtained by running cdd or other program and 

6 The Sub-Cone method 

We give here a method that was applied to only one cone, i.e. the hypermetric cone on seven 
vertices HYP-j. The cone HYPj is defined by 3773 facets called hypermetric facets which form 
14 orbits under the group Sym{7) (see [B70] and [Ba99]) 



1. 



We have computed the facets of the cone CUT^. There was 36 orbits of facets [Gr90]. 10 
of them are hypermetric. 



2. The cone HYP-j is decribed by 14 orbits of hypermetric facets, 10 already appear in the 
hst of 36 orbits of facets 

3. We define 26 subcones 

Hp = {xe HYPj and p{x) < 0} (1) 

for every representative p of the 26 orbits of non-hypermetric facets. We remove the 
redundant facets by hnear programming (see [Fu95]) 

4. The non-redundant description of Hp are simphcial cones which were easily computed, so 
giving us description of extreme rays of HYP-j. 

The extreme rays of HYP^ belong to two classes: the cuts coming from CUT^ and the 26 found 
extreme rays. Theses extreme rays are adjacent only to cuts, so the adjacency decomposition 
method is not useful here: The skeleton graph of HYP-j without the cuts is totally disconnected. 

The success of the method rely on several factor: the knowledge of facets of CUT-j which is 
not something easy [Gr90] and the fact that the cones Hp are disjoint making the computing 
process reasonable. The cuts form a complete graph in HYP^ but this fact solely does not 
explain why the Hp are disjoint since MET-j which is in the same situation as HYP-j does not 
share this property. All attempt to use this method in other contexts were inefficient. 
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